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We introduce a new quantum transport formalism based on a map of a real 3-dimensional lead- 
conductor-lead system into an effective 1- dimensional system. The resulting effective ID theory is 
an in principle exact formalism to calculate the conductance. Besides being more efficient than the 
principal layers approach, it naturally leads to a 5-partitioned workbench (instead of 3) where each 
part of the device (the true central device, the ballistic and the non-ballistic leads) is explicitely 
treated, allowing better physical insight into the contact resistance mechanisms. Independently , we 
derive a generalized Fisher-Lee formula and a generalized Meir- Wingreen formula for the correlated 
and uncorrelated conductance and current of the system where the initial restrictions to ballistic 
leads are generalized to the case of resistive contacts. We present an application to graphene 
nanoribbons. 

PACS numbers: 72.10.Bg, 73.23.-b, 73.63.-b, 73.40.Cg 



I. INTRODUCTION 

The problem to describe nanoscale electronic 
transport^ from first principles still remains a 
formidable challenge. Although powerful formalisms 
have been developed in the last years^^ the cornerstone 
is still represented by an efficient coupling between 
electronic structure theories with appropriate modelling 
of the quantum transport problemi&i In the framework 
of the Kubo formalism^ 2 - very efficient methods, such as 
the MKRT approach^ have allowed a wealth of appli- 
cations even to technological systems.— The Landaucr 
formalism and its NEGF extension-^ recently witnessed 
more intense developments which have led to the setup 
of a standard model, often referred as principal layers 
approach^ In this approach, the workbench model is 
assumed to be a 3-partitioned system^ constituted by 
a central region sandwiched between two semi-infinite 
leads. The latter, assumed ballistic and at partial 
equilibrium, only inject and harvest electrons into the 
central region where all the processes affecting the 
conductance are assumed to take place. Those typically 
include contact resistance, scattering by impurities 
and defects, incoherent transport electron-electron 
and electron-phonon scatterings. Once provided the 
electronic structure of the system by an appropriate 
(ab initio or semi-empirical) theory, the consecutive 
quantum transport problem can be solved by e.g. the 
calculation of the Green's function of the central region 
in presence of the effect of the leads, represented by 
the self-energies of semi-infinite periodic systems. The 
conductance can then be calculated via the Fisher-Lee 
or the Meir- Wingreen formulas^ involving the Green's 
function and the leads injection rates. 

One may notice that in this standard model the sepa- 
ration between the central region and the ballistic leads 
appears somehow arbitrary and unphysical. Indeed, in 



order to correctly describe the contact resistance, the 
central region should contain not only the conductor un- 
der study (e.g. a molecule or a nanodevice), but also 
some layers in reality belonging to the leads. Conver- 
gence should be checked by increasing the central region 
size and thus the number of states in the problem. Its 
computational resolution is heavier since it deals with a 
number of channels much greater than the true channels 
of the central device. 

In this work we introduce the notion of effective chan- 
nels as the states through which the current flows up 
to the central device. The number of these channels is 
upper-bounded by the number of states of the central 
bottleneck. All the leads' states orthogonal to the states 
of the effective channels do not effectively participate to 
the conductance and can be safely disregarded, acting 
as a prcfiltering of the initial problem. This is a consid- 
erable simplification with respect to the principal layers, 
where the number of the channels is usually much greater 
than the number of channels in the device. The numer- 
ically efficiency is consequently much improved. More- 
over, the effective channels can be viewed as associated 
to an effective 1- dimensional system into which the real 
3-dimcnsional physical system is mapped (Fig. [I}. This is 
a way to restore the natural dimensionality of the quan- 
tum transport problem, which is truly ID. The resulting 
ID effective theory is an in principle exact formalism to 
calculate the conductance. The main advantages of the 
ID approach are: 

• It is formally exact, and its numerical implementa- 
tion requires no approximation. 

• It reduces the size of the numerical problem. 

• It is particularly efficient for nanocontacts. 

• The calculation scales as N tot x iV 2 hanncls , with N tot 
the full size of the initial Hamiltonian and A^hanneis 
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the number of effective channels. 

• ft is independent from the level of theory used to 
calculate the electronic structure, whether with or 
without correlation. It relies on the preliminary 
knowledge of the Hamiltonian and eventual inter- 
action self-energies. 

• It is independent from the type of basis set used in 
the previous electronic structure calculation. 

Indeed, the formalism only requires a complete -/V c h a nneis 
set of independent states localized at the bottleneck, 
as well as how H acts on states, as implemented in com- 
mon DFT codes. Afterwards, there is absolutely no re- 
quirement on the basis set that can be atomic orbitals, 
Wannicr functions, non-orthogonal or non-maximally lo- 
calized, even plane waves. In fact, another important 
result of this formalism is that it proposes its own recur- 
sion states basis as the best suited to solve the quantum 
transport problem. 

Present limitations and drawbacks of the approach are: 

• It has one step more (the recursion) than a princi- 
pal layers calculation. 

• Although a generalization to the feii -dependent 
transport is possible, the approach is inefficient on 
non-bottleneck, such as planar, geometries. 

• For transport at finite voltage, the effective chan- 
nels must be recalculated at each different bias. 

Independently, in this work we will also introduce a 
physically more intuitive 5 -'partitioned (instead of 3-) new 
quantum transport workbench model (Fig.s[2]and|3]). This 
is composed by the true central conductor device, the left 
and right sections of non ballistic leads — which con- 
tain and isolate contact resistance mechanisms — and 
finally the ballistic semi-infinite leads. We will derive the 
exact uncorrelated and correlated conductance formulas 
associated to the general workbench where a section of 
the leads is non-ballistic. For the uncorrelated case, we 
will derive a generalized Fisher-Lee formula to be asso- 
ciated with the 5-partitioned workbench. For the case 
of correlated transport we will derive a generalization of 
the Meir-Wingreen formula to resistive nanocontacts. In 
both cases the effect of the contact resistance is exactly 
taken into account and contained into renormalized lead 
injection rates T and scattering functions S <> . The new 
workbench allows a valuable insight on the origin of re- 
sistance into the system and provides a clear analysis of 
the different resistance mechanisms. While in the princi- 
pal layers all the resistance is contained in the extended 
central region Green's function Gc, here the contact re- 
sistance can be directly read from T and E <> whereas 
other more internal mechanisms — i.e. density of states 
effects, scattering by impurities or defects, and e-e and 
e-ph scattering effects — can be read from the Green's 
function G c of the true central part. The application of 
the generalized Fisher-Lee and Mcir-Wingrccn formulas 
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FIG. 1: Mapping of a real 3D device (top: hydrogen molecule 
in between gold leads) into the effective ID system (bottom: 
effective atomic chain). The effective channels arise from the 
central device (here the hydrogen molecule) and pursue into a 
non-ballistic section (blue, violet and red pseudoatoms) , until 
they achieve an asymptotic ballistic behaviour (yellow). 



is not at all restricted to an effective ID problem and 
can be straightforwardly implemented also in a ordinary 
3D quantum transport geometry (Fig. [3]). Any actual 
principal layer code can be modified in order to refer to 
the 5-partitioncd workbench and the associated quanti- 
ties. With respect to the principal layers, the 5-partioned 
workbench model seems to present only advantages: 

• It allows a clear analysis of the resistance by disen- 
tangling contact resistance from other mechanisms. 

• The initial problem is divided into several indepen- 
dent parts, with reduced calculation cost. 

• It is very efficient on the effective ID system. 

We will show an application to graphene nanoribbons 
in the tight-binding approach. 



II. DEFINITION OF EFFECTIVE CHANNELS 

Let's first provide a heuristic approach to the notion of 
effective channel. A generic nanoscale quantum transport 
system is always characterized by a bottleneck of a few 
N a atoms, each contributing with say max A^ Q orbitals. 
For example the bottleneck is represented by a hydrogen 
atom in Fig. [I] a nitrogen atom in Fig. [3] and by 6 carbon 
atoms in Fig. SJ The conductance of the whole system 
can never exceed -/V c hanncis = N a N a quantum of conduc- 
tance. In the ballistic case the conductance has the char- 
acteristic steplike integer profile (Fig. [4] dashed line) , at 
most summiting to N a N Q . In the other cases, the con- 
tact resistance always reduces the conductance (Fig. 21 
continuous line). The bottleneck represents an upper- 
most bound. So we can say that there are at most N a N a 
channels effectively contributing to the conductance, the 
rest being idle. This notion of effective channels can thus 
be exploited to simplify the problem. 
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Let's now find a general condition for a state to con- 
tribute to transport to be used as definition of effective 
channel. We consider a quantum transport system sepa- 
rated into a true central device c (for example "only" the 
hydrogen molecule of Fig. [IJ coupled to left I and right 
r contacts. The Hamiltonian can be written: 

H = H c + Hi + H r + H tc + H cl + H cr + H rc , (1) 

where H c , Hi and H r are the Hamiltonians of the cen- 
tral device, the left and right contacts respectively, and 
H c i = H( c and H rc = H\ r the coupling of the central de- 
vice to the left and right contacts. The absence of direct 
coupling between the two contacts is here a fundamental 
requirement. An electron must pass across the central 
device to go from one contact to the other. We name 
S c , Si and S r the Hilbert spaces of states of the central 
device, left and right contacts. 

We define the effective channels space of lead t (t = I 
or r) by 

Sf = Span{iT t "#tc|<M} V0 C g S c , Vn 6 N. 

Sf C St is a Hilbert subspace of St containing only 
the states coupled to the central device. On the other 
hand, its orthogonal complement <S t (such that St = 
S^ © Sf- , direct sum) contains only states which are H- 
disconnected both from S c and Sf. The application of 
H to a state of Sf- still belongs to Sf- . This means that 
an electron from the reservoir t but in a state belonging 
to Sf~, will stay and evolve in Sf- without contributing 
to a current across the central device. The calculation of 
the conductance is not affected by Sf- that can be safely 
neglected. Thus the conductance of the real system is 
that of a simpler effective system living into S e ^, defined 

S eB =Sf +S c + Sf. (2) 

III. ID EFFECTIVE SYSTEM 

The next step is to demonstrate the existence of this ef- 
fective system by a correct choice of its basis. In practice 
we will provide a Gram-Schmidt or, equivalently, a Hay- 
dock recursion^ algorithm to build an orthonormal basis 
set {il> n } for S e $. The representation H e ^ of the original 
Hamiltonian on {ip n \ will result at the same time. 

To start with, let's restrict to the simplest (scalar) case 
as in Fig. [1] where the central bottleneck is a single atom 
with a single orbital, say 4> c . Let's first build the right Sf 
effective channel space and its basis set {t/jn}- The first 
element ipi of the basis is given by bi\ipi) = H rc \4>c)- b\ is 
chosen as normalization factor for ipi . Next we calculate 
a i = ( , 4>i\Hr\'4>i) ■ We then calculate the second element 
by 62^2) = H r \tpi) — OilV'i)- "02 is orthogonal to 
and normalized by 62. At the next and all the following 
steps we iterate the same procedure, a n = (ip n \H r \ip n ) 
and b n +i\ip n +i) = H r \ip n ) - a n \ip n ) ~ Kli'n-i}- This is 
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FIG. 2: The conductance can be calculated by the traditional 
Fisher-Lee formula and 3-partitioned workbench, with ballis- 
tic leads and extended (molecule + non-ballistic leads) central 
region (scheme above); or the generalized Fisher- Lee formula 
and the 5-partitioned workbench, with central device, non 
ballistic, and finally ballistic sections of the effective channels 
(scheme below). 

an implementation of the standard recursion method^ 

H r \ip n ) = a n \ip n ) + b*\ip n -i) + b n+ i\ip n+ i). (3) 

In conclusion we end with an orthonormal basis {ipn} for 
Sf. With increasing n, the state ip n is a linear com- 
bination of real orbitals belonging to atoms deeper and 
deeper in the contact ^ The recursion can be stopped at 
an n = N where the coefficients a n , b n saturate and con- 
verge to an asymptotic regime, et^, b^ (see Fig. [5]). From 
this point the leads are consequently ballistic and associ- 
ated to states achieving a maximum spread into the con- 
tact region. Notice that this practically and numerically 
recovers the Landauer's paradigm of reservoir^ initially 
formulated as an ad hoc hypothesis. On this basis set the 
Hamiltonian is tridiagonal, with onsitc = a n and 
(only) first neighbours hopping coefficients H ( f n _ 1 = b n . 
It can be seen as associated to an effective ID pseu- 
doatomic chain (Fig. [IJ. H e $ is in fact the same origi- 
nal Hamiltonian but represented on an orthonormal basis 
where it is tridiagonal. Notice that the algorithm does 
not require to store the basis elements tp n , but just only 
the a n and b n coefficients on the basis (see Appendix 

EH) . 

The procedure can now be repeated for the left I con- 
tact ({ipn}, n < 0). And it can be generalized to the 
matrix N a N a ^ 1 case (see Appendix I A 2[) . In this case 
a n and b n are replaced by block matrices A n and B n of 
size N a N a . The Hamiltonian still owns a ID structure 
but it is now tridiagonal by blocks in the recursion basis: 

N N 

H^=H C + J2 A "+ E ( B n + Bi). (4) 

n=— N n=-N 



IV. GENERALIZED FISHER-LEE FORMULA 

We now focus on the calculation of the conductance 
of the effective ID system. The principal layers ap- 
proach would reorganize the system into 3 new regions 
(Fig. [U top) : a left ballistic lead region L where a n sat- 
urates to the constant asymptotic value a_oo; a right 
R ballistic lead region a +oa \ and an extended central 
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into the Fisher-Lee Eq. ([5]) and using some fundamental 
projector relations^ we get at (see Appendix IB 4[) 



FIG. 3: Illustration of the 5-partitioned workbench for a 
generic 3D quantum transport system, with central device 
c, non ballistic A, p, and finally ballistic sections L and R of 
the leads. F ; / r are the contact resistance renormalized lead 
injection rates, and g p /\ the Green's functions of the leads 
non-ballistic sections A and p calculated as they were dis- 
connected from the central region c and in presence of the 
external ballistic leads L and R. 



region C containing the true central device c and also 
the two non-ballistic sections of the leads, n = 1, . . . , N 
and n = — 1, . . . , — N. Referring to this standard work- 
bench, we can calculate the Green's function g^ and g R 
of the semi-infinite periodic leads, the associated self- 
energies S L = H C l9lHlc, = H CR g R H LR and the 
injection rates Tl, T r . We can then calculate the re- 
tarded/advanced Green's function Gp° of C in presence 
of the leads L and R by G c (z) = (z - H c - £l - ^r)' 1 - 
The conductance is then calculated by the standard 
Fisher-Lee formula, 



C(z) 



2e 2 



Tv[T L {z)G c {z)T R {z)G c {z)) 



(5) 



The drawback of this procedure is twofold. On one hand 
there is no way to analyze separately the role of the true 
central device and of the non ballistic part of the leads. 
On the other hand the size of the matrix G(z) can in- 
crease rapidly and the matrix inversion can become dif- 
ficult or impossible, even for a tridiagonal hamiltonian. 
This restricts the applicability of the method to systems 
where the resistance is localized at the very vicinity of 
the nanocontacts. 

In this work we instead propose to keep the original 
natural separation into true central region c and the leads 
I and r. Into I and r we identify the ballistic regions 
L and R and the non-ballistic sections A and p (Fig. [5] 
bottom and Fig. [3|. By using the projector P c on c, we 
single out the Green's function of the true central device 
G c = P c GcP c - We then define the Green's function g\ = 
[z — P\(Hc + ^l)P\}^ 1 (a similar expression for g p ). This 
is not the projection of the propagator into the section A, 
9\ 7^ G\ = P\GcP\- Instead it is the propagator in the 
section A calculated as if A were disconnected from the 
center c but connected to L. Injecting these definitions 



C(z) 



2e 2 



where 



and 



t r (z) 



H cl g a x (z)T L (z)gi(z)H lc , 
H cr ~g r p (z)T R (z)g a p (z)H rc , 



(6) 



(7) 



(8) 



G c = {z-Hc-Ex-Hp)- 1 , 

= H c ig x (z)Hi c , 
S p = H cr g p (z)H rc . (9) 

Formula ((6]) has a Fisher-Lee like form but involves dif- 
ferent quantities (for instance T;/ r ^ i[S£, — £^ p ]). 
It now refers to a workbench where the Green's function 
Gc of the extended central region is replaced by the more 
significative Green's function G c of the true device un- 
der study. The injection rates T L / R of ballistic leads are 
replaced by contact resistance dressed renormalized injec- 
tion rates Tu r which refer to both the ballistic L/R and 
the non ballistic X/p sections of the leads. In the princi- 
pal layers approach all resistance mechanisms are local- 
ized within the extended central region C and considered 
in Gc- Here contact resistance is separated from other 
mechanisms, localized in the non-ballistic sections A and 
p and transferred into T where it is taken into account 
via g. The contact resistance can be read directly from T. 
As we will see in the example, if T(E) — at a given E, 
this will provide conductance whether or not there is at 
E an available channel in the central device. Therefore 
the generalized Fisher-Lee formula allows a more clear 
interpretation of resistance mechanisms. Notice that the 
Ti/ r depend only on the electronic structure of the con- 
tact and on its coupling to the central device. Thanks 
to recurrence relations the calculation of the ID g and 
hence of T can be carried out in a very efficient numerical 
way (see Appendix IB 4[) wrt the principal layers approach 
where the calculation of Gc can be cumbersome. More- 
over, the generalized Fisher-Lee formula is not restricted 
to the ID-effective theory presented here. It can be ap- 
plied also to real 3D systems, in place of the ordinary 
Fishcr-Lcc, provided the extended central device is split 
into two non-ballistic leads regions and the true central 
device (Fig. One might think for example about a 
molecular junction where our formula could be used to 
project the whole transport problem into molecular or- 
bit als. 



V. EXTENSIONS TO NEGF AND 
CORRELATED TRANSPORT 

The 5-partitioncd workbench is particularly convenient 
in the case of correlated transport within NEGF. Start- 
ing from the Meir-Wingreen formula^ for the current, 
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FIG. 4: Conductance of a metallic zigzag <5.0> graphene 
nanoribbon (left) of length L — 2.1 nm connected to graphene 
2D seminfinite sheets. The zero of the energy corresponds to 
the charge neutrality point of graphene i.e. to the Dirac point. 
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FIG. 5: Convergence of the 6 (i = 1, . . . , 6) a l „ eigenvalues of 
the matrix A as function of the iteration Ti- 



the non-coherent term can be separated from the coher- 
ent. For the coherent term we end again to the gener- 
alized Fisher-Lee Eq. ©. For the non-coherent current 
we can derive a generalized Meir-Wingreen expression of 
the form: 

*r coh = I Tr[E<G^S> „G a - E>G£E< , (10) 

where S< > r are the in/out scattering functions related 
only to correlations (e-e or e-ph), while E^ are the 
in/out contact resistance renormalized lead t scattering 
functions, T, t (z) = H ct g^(z)T, T {z)g r T {z)H t c, w ith T , T (z) 
calculated in the ballistic region (see Appendix IB 6|) . No- 
tice that E t 7^ E T , with E T defined in Eq. (j9]). The equi- 
librium relations E< = /f D f 4 and E> = (1 - jf D )f t 
(as well as T t = [Ef — Ef ]) still hold also for the renor- 
malized quantities. Contact resistance is now physically 
separated into E t with respect to the resistance raising 
from e-e and e-ph scattering mechanisms, associated and 
localized into the true central device and G c . This more 
faithfully represents the workbench ideal assumption of 
lost-of-coherence effects only within the true central de- 
vice, with leads assumed as everywhere perfectly coher- 
ent. 



VI. APPLICATION TO A GRAPHENE 
NANODEVICE 

The method has been implemented into a computer 
code^i and was applied to graphene nanoribbons coupled 
to graphene sheet s 20 ' 23 using a tight-binding electronic 
structure. As shown in Fig. |4] and in Ref.— the conduc- 
tance exhibits Fabry-Perot oscillations shorting even to 
some 0s close to E = 0. Although the nanoribbon is 
metallic (A C (E ~ 0) ^ 0) the transmission of electron 
waves across the contact is blocked at these energies. It 
can be shown that the renormalized T t (E) is zero at those 
points, whereas Tt keeps finite. Thank to the new for- 
malism, our analysis shows that the conductance is 
a pure effect of contact resistance and it can be inter- 
preted as a diffraction effect at the contact constriction 



(see ReL2°-). Notice the high accuracy of the calculation: 
the maxima of the Fabry-Pcrot oscillations, in principle 
exactly equal to 1, are found equal tol±10~ 3 ~10~ 4 . 
Furthermore, at difference with respect to the principal 
layers approach as implemented e.g. in WanT— or other 
codes, the conductance as function of the energy does 
not present spurious structures and spikes shorting to 0. 
Notice also the straightforward, well defined and reliable 
convergence criterion. In Fig. [5] we show the typical con- 
vergent behaviour of the a n coefficients as function of 
the iteration n (the 6 eigenvalues of the A n matrix in the 
present case). The code automatically checks the con- 
vergence and stops the recursion when \a n — a n -i \ < SE 
and then calculates the conductance. To get the 10 -3 eV 
accuracy we needed n = 600 iterations. The case studied 
is however one of the most unfavourable, due to a diver- 
gence of the electron wavelength at E = 0. We expect a 
much quicker convergence for less critical systems, where 
contact resistance is localized closer to the nanocontacts. 



VII. CONCLUSIONS 

We have introduced a formalism for quantum trans- 
port at resistive nanoscale contacts which relies on the 
introduction of an effective ID system and its associated 
5-partitioncd workbench model. Moreover, we generalize 
the Fisher-Lee and the Meir-Wingreen formulas to the 
case of non-ballistic contacts. The formalism is physi- 
cally more intuitive and numerically more efficient than 
the standard principal layer approach. By its versatil- 
ity and its theoretical generality, our formalism allows 
to handle unexplored quantum transport problems, like 
the exact influence of the experimental contact geometry. 
From a more fundamental point of view, our work offers 
new perspectives to discuss the existence of Landauer's 
reservoirs for real systems. 
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Appendix A: Mapping into the effective ID system 
and matrix recursion 



1. Construction of the block spaces 

Let us use the convention where the recursion into the 
right terminal t — r is indicated by positive recursion 
indices n > 0, while for the left contact t = I by negative 



indices n < 0. Let us define the space S 1 



"11 



Sf=Span{H rc \(P c )}, 



(a similar expression holds for Sff = Span{J?i c |0*)}). 
The application of H r to generates states that can be 
decomposed on itself and on an orthogonal space that 
we name . Successive application of H r to states of 
generates states that can be decomposed on S^, S^f 
and a new orthogonal space 1S3 . In general, application 
of H r to Sf-p generates states that can be decomposed on 
S^_ x , Sfp and a new orthogonal space S^_ x . Indeed the 
Hermicity of H r prevents the coupling between spaces 
that are not successive order. 

The effective spaces for the left and right terminals are 
hence given by 



S, 



eff . 



n>0 



eff 



ri<0 



eff 



Let us now define the orthogonal projector P n = P£ = 
P% on the space S%&. A„ = P n HP n is the restriction of 
the Hamiltonian H to the space and B n = P n HP n _i 
represents the coupling from S^~_ 1 to S%f. The coupling 
from to S^_ x is given by B\. One can see that 



H r = *£A n + J2(Bn + Bt) (Al) 

n>0 n>l 

H i = E A » + E ( B n + Bi) (A2) 



n<0 n<— 1 



and 



H rc + H cr — B\ + B[ 
Hie + H c i = B-i+B_ x 



(A3) 
(A4) 



Therefore the Hamiltonian restricted to the effective 
space S e % = Sf& + S c + can be written: 

N N 

H e ff=H c + A -+ E (B n + Bl). (A5) 



2. Matricial recursion 

In the general case the subspace spanned by H tc \4>c) 
with <p c any orbital in the central device, is of dimension 
N a N a , that is the characteristic dimension associated to 
the smallest bottleneck in the central device. In princi- 
ple, we can start the recursion from the bottleneck itself 
by finding N a N Q linearly independent states $ c (i) with 
i = I, . . . , N a N . In any case, the dimension of S± or 
«Sjf , as well as all successive spaces S%f?, is strictly upper- 
bounded by N a N . The application of the first iteration 
recursion step leads to N a N„ linearly independent states, 
say ^i(i) with i = 1, . . . ,N a N Q . The same for all suc- 
cessive iterations, ^ n (i). Then a„ and b„ are replaced 
by matrices A n and B n of dimension N a N J£ The re- 
currence relations have an analogous form. We consider 
here the case for the right terminal r and the positive 
recursion n > 0. The projector on the subspace S!fp can 
be written as 



p„=x;i*„(*)><*»(i)i 

i 

Then one get from Eq. (|A1[) : 

^r|*n(0)=X)[A»(j.i)l*nC/)> 



(A6) 



with 



Bi(t,j)\^ n ^(j)) + B n+1 (j,t)\^ n+1 (j))} (A7) 



(*n(i)l*m0')> = <WAj 

except for n = |*o(*)> = (A8) 



and 



A n (i,j) = (* n (i)\H\* n (j)) (A9) 
B n+1 (i,j) = (* n+1 (i)\H\* n (j)) (AIO) 

The procedure to compute the matrices A n and B n is 
quite similar to that of the scalar case. At step n one 
knows A m and B m+ i for all m < n. One stores also the 
components of \^ n ) and |\I/ n _i} in the real space basis. 

Then A n (i,j) is calculated from equation (| A9|) . One 
then compute the components in the real space basis of 
K+i(i)>: 

K + iW} = E B "+ i ^ i )i*™+iW) = 

3 

H r \*n(i)} -EKO'.OhMj)) - 
3 

St(i,i)|*„_ 1 (j))] (AH) 

with the overlap matrix 

Sn+l&j) = <*«+l(*)l*n+l0')> = (Bi +1 B n+1 )(l,j) 

(A12) 
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Once the overlap matrix S n +i(i,j) is calculated, we can 
construct the orthonormal basis of the ^ n+ i(j) states 
whose components on the basis of the ^ n+ i(j) states are 
given, as well as of course on the real space basis. Once 
the <i>„ + i(j) states are calculated, B n+ \{i,j) is calculated 
from Eq. (|A10[) . We then pass to the next step of recur- 
sion. 

There is a freedom on the choice of the orthonormal ba- 
sis of the subspaces S^_ 1 which is generated by \I/ n +i(j) 
for all j, or equivalently by ^ n+1 (j). Indeed the matrices 
A n and BnB^ are defined up to a unitary transformation 
which depends on the precise choice of the vectors \& n (i). 
This unitary transformation is equivalent to the freedom 
in the choice of the phase of the b n coefficients in the 
scalar case. 

The most time consuming part is the computation 
of the Hamiltonian of the effective ID system. This 
amounts essentially to compute 2N a N scalar recursion 
procedures, where N a N a is the dimension of the spaces 
S^P corresponding to the size of the blocks of the tridi- 
agonal Hamiltonian H e $. Since the only operations done 
on the real space basis are scalar products, the whole pro- 
cedure exactly scales linearly with the size of the initial 
Hilbcrt subspace. That is, the proposed algorithm is an 
O(N) method. Once the Hamiltonian Eq. (|A1[) is com- 
puted, the calculation of the operators Tt(z) and S t (z) 
is relatively fast. 



Appendix B: Derivation of the generalized 
Fisher-Lee and Meir-Wingreen formulas 

In this appendix we present two different schemes 
of derivation for the generalized Fisher-Lee and Meir- 
Wingreen formulas. The first scheme is an ordinary 
matrix derivation, while the second uses projector tech- 
niques. 



1. From the 3- to the 5-partitioned workbench 



The original Hamiltonian Eq.([T]) 

Hi Hie 
H = I H c i H c H cr 
H rc H r 



(Bl) 



is written in the physically intuitive 1— c— r 3-partition of 
the system, with c the true central region, and I and r the 
external leads containing both ballistic and non-ballistic 
sections 

In the principal layers approach, the Hamiltonian is 
divided differently and is rewritten as 

Hl Hlc 
I H CL H c Hcu. 
Hrc Hr 



that is using a different L-C-R 3-partitioned scheme 
where the system is divided into left and right ballis- 
tic leads L and R, and a central extended region C that 
now contains the non ballistic sections of the leads. The 
associated Green's function is 





t G L 


Glc 


Glr 


G = 


Gcl 


G c 


Gcr 




\ Grl 


Grc 


Gr 



and we use the ordinary principal layers formulas to cal- 
culate the conductance (T = L or R): 



9t = {z- H t ) 1 

Tit = HctQtHtc 

Gc = {z-H c -^l- Sh)- 1 



(B2) 



2e 2 



C = —Tr[T L G c T R G c ] 



Starting from the principal layers scheme, we now di- 
vide C in 3 subsections: the left non-ballistic lead section 
A, a true central region c and a right non-ballistic lead 
section p (see Fig. [2] and Referring back to the 1-c-r 
scheme, 1 is the left lead comprising both the ballistic L 
section and the non ballistic section A, and r is the right 
lead comprising both the non ballistic section p and the 
ballistic section R. We end it up with a 5-partitioned sys- 
tem L-A-c-p-R (see Fig. [3]) and the associated Hamilto- 
nian 



( HL 

H\l 



H = 



\ 



Hlx 
H x 
H c \ 









H\ c 

H r 






H r , 



\ 



H nr H n 



H 



Rp 







HpR 

Hr J 



ED we 



Comparing with the original Hamiltonian Eq. 
see for example that 

Hd = ( H cX ) , 



and so on, so that in the next we will confuse H ct with 
H CT (t — I or r and T = Aorp). whenever it is clear from 
the context in which space we are working. In principle 
we start from the principal layers approch and we will be 
always working in the space C. 

The extended central Green's function Gc is rewritten 
in its X—c-p components: 



G c 



G\ G\ c G\ p 
G c \ G r _ G 
G 



G 



'■/> 



'PC 



G B 



The injection rates Tl and Tr, which live in the C space, 
have non-zero elements only in the regions A and p re- 
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spcctivcly, and are of the form 



r A o o 
ooo 
ooo 

ooo 
ooo 
o o r„ 



(B3) 
(B4) 



like also E^ and E#. 



2. Generalized Fisher-Lee formula 

Starting from the ordinary Fisher-Lee formula refer- 
ring to the principal layers approach and the L-C-R 
workbench, 

2e 2 2e 2 
C=— T=— Tr{r L G r c r R G a c ], 



(B5) 



using Eq.s (|B3j) and (|B4j) we rewrite T as 
T = Tv[T x G r Xp r p G a pX ]. 



G\ p and G p \ arc determined cxpliciting Eq. (|B2[) on the 
A-c-p scheme, 



H\ + Ea #a, 
#cA #c 
i? 






#cp 



From 



(z — H\ — Ea)Gac — H\ C G C = 
—Hp C Gc + (z — Hp — Y>p)Gpc = 
—H C \G\ C + (z — H C )G C — H C pGpc 



= 1. 



we get 



Gac = gxHx c Gc 

G pc — QpHpcGc 

G c = {z-H c -t x -±p)-\ 
where we have defined 



,9a = {z-H\- 

9 P = {z-Hp- 

Ea = H c \g\H\ c 

Ep H C pQ pH p C . 



E„) 



(B6) 



Notice that <? r ^ G T (where r = A or p). The <? T can 
be physically interpreted as the Green's function of the 
non-ballistic lead section r (see Fig. \S§ as disconnected 
from the central region c (H CT = H TC = 0) but in contact 



to the ballistic lead region T and under the effect of its 
E T (E t ). 
Now from 

(z-H x - Ea)Ga p - HxcGcp = 
—Hp C Gcx + (z — H p — E p )GpA = 
(z — Hx — Ea)Ga — HxcGcx = 1 
—Hp C G C p + (z — Hp — Yip)Gp = 1, 



we get 

Gap = 9xHx c G C p 

GpA = 9pH pc G c x 

Ga = <?a + gxHxcGcx 

G p = g p + 9pH pc G C p. 

And finally from 

—HcxGx + {z — H c )G c x — H cp Gpx 
—H c \G\ p + (z — H c )G cp 



(B7) 
(B8) 



H C pGp 



- 
0. 



we get 



GV 



G c H c xgx 



Gcp G c H C pQp, 



which replaced into Eq.s (|B7[) and (|B8|) . provide 

Gap = gxHxcG c H cp gp 
GpA = g P Hp C G c H c xgx, 



that we can finally insert into Eq. (|B5I) to get the trans- 
mittancc 

T = Tr\Txg\Hx c GlH C pg r p T pg^HpcG^.Hcxg'i] 
= Tr[H c xg^*9xH*cGlH cp g p r p g p H pc Gc]- 

Defining the contact resistance dressed renormalizcd lead 
injection rates V, 



T; — Hcxg^T LfjxHxc 
T r = H cp g p T R g a p H pc 



(B9) 
(B10) 



we get the final result for the generalized Fisher-Lee for- 
mula 

r = Tr[f/G^f r G°]. 
Finally we notice that 



as it can be checked from the definitions Eqs. (|B10[) and 
(|B6|) . We can however define a new E t ^ E r , 



E t = H CT g°Z T g r T H T 
such that the relation 

f t = i[E>-E t °], 

holds. 
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3. Fundamental relations for the resolvant 

We now present the projector methodology developped 
by Zwanzig and Morii^ that we will use later to provide 
a different derivation of both the generalized Fisher-Lee 
and Mcir-Wingrecn formulas. We will first rederive some 
fundamental relations that we will use in the next sec- 
tions. Let's consider a space V and its associated projec- 
tor P. Let's define also its complementary Q, that is the 
projector associated to the complementary space Q. By 
their definition, the following relations hold for P and Q: 

P + Q = l, 
P 2 = P, 

Q 2 = Q, 

PQ = QP = 0. 

We now consider a not necessarily Hcrmitian operator 
H (which can even depend on z although we will drop in 
the following this dependence) and its associated Green's 
function G(z) = (z — H) . Starting from the identity 

{z-H)G= 1, 

and multiplying from the left by P and from the right by 
Q, we get the relation 

P(z - H)GQ = PQ = 0. 

We now insert the identity operator (P + Q), 

= P(z - H){P + Q)GQ 

= P(z - H)PGQ + P(z - H)QGQ 
= P(z - H)PPGQ - PHQGQ, 

and we finally arrive to the fundamental projector rela- 
tion 

PGQ = [P{z - H)P]~ 1 PHQGQ. (Bll) 
Starting from the identity 

QG{z - H)P = 0, 
we arrive to the other Zwanzig-Mori projector relation 
QGP = QGQHP[P(z - H)P]-\ (B12) 

4. Generalized Fisher-Lee formula by the projector 
scheme 

With respect to the X-c-p regions scheme, we define 
the projectors P\, P c and P p which undergo the following 
relations 

P x +P c + P p = l, (B13) 
Qx = Pc + P P , (B14) 
Q P = Px + Pc, (B15) 
QxP P = P P , (B16) 
QpPx = Px- (B17) 



For a tridiagonal Hamiltonian H' c (here we use the 
convention to include into H' c not only the Hamiltonian 
He associated to the region C, but also the lead L and 
R self-energies, H' c = He + £l + £.r), the only non-zero 
crossed terms are 

Hxc = PxH' c P c , (B18) 
H cX = P c H' c P Xl (B19) 
H cp = PcH'cPp, (B20) 
H pc = P P H C P C , (B21) 

while the only non-zero T elements are: 

Tl = PxTlPx = r Aj (B22) 
r R = P p T R P p =T p . (B23) 

We now start from the ordinary Fishcr-Lce formula 
referring to the principal layers approach and the L-C-R 
workbench, 

2e 2 2e 2 
C=—T=—Tr[T L G r c T R G a c }. 

We rewrite the transmittance T using Eqs. (|B22[) and 
(IB231) 

T = Tx[P x V L P x G r c P p T R P p G c ]. 

Using the cycling property of the trace, 

T = Tv[T L P x G c P p T R P p G a c P x ], 

and the projector properties Eqs. (|B16j) and (|B17j) . wc 
arrive to 

T = TiiT L {P x G r c Q x }P p T R {P p G a c Q p }P x }. 

Now we use the Zwanzwig-Mori relation Eq. (|B11[) 

T = Ti\r L {[P x {z-H' c )P x )- x P x H' a Q x G r c Q x }P p T R 
{[P p (z - H' c )P p ]- l P p H' c Q p G a c Q p ) P X ]. (B24) 

We now define the Green's functions g, 

g x = [P x (z-H' c )P x ]- 1 , (B25) 
g p = [P p (z-H' c )P p ]- 1 . (B26) 

Notice that g x ^ P x (z - H' C )~ 1 P X and g p ^ P p (z - 
H' c )~ 1 P p . The g x (g p ) can be physically interpreted as 
the Green's function of the non-ballistic lead section A 
(p) (see Fig. [3]) as disconnected from the central region 
c but in contact to the ballistic lead region L (R) and 
under the effect of its (Si?). Replacing the definition 
Eqs. (|B25l) and ([B26) into Eq. ([B24| . we get at 

T = Tr[T L f x P x H' c Q x G r c Q x P p T R ~g a p PpH' c Q p G a c Q p P x ]. 

Developping Q x and Q p according to Eqs. (|B14[) and 
(|BT5|) and using Eqs. ([B18)l and (|B%|) for H Xc and H pc , 
we get 

T = Tr[T L glH Xc P c G r c P p r R g a p H pc P c G a c P x ]. 
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Using the projector relations P c = P C Q P and P c = P C Q\, 

T = Ti[T L ~g r x H Xc P c {Q p G r c P p }T R g a p H pc P c {Q x G a c Px}l 

and the Zwanzwig-Mori relation Eq. (|B12[) . we get 

T = TT[r L ~g r x H Xc P c {Q p G r c Q p H' c P p [P p (z - H^P,}- 1 } 
T R ~g a p H pc P c {QxG a c Q x H' c P x [P x (z - ^)P A ]~ 1 } ]. 

Using again the definition the Green's functions g and 
the project relations, 

T = Tv[T L g r x H Xc P c G r c P c H' c P p ~g r p 
T R ~g a p H pc P c G a c P c H' c P x g a x }, 



and the definitions of H cp and H cX Eqs. (|B19|) and (|B20I) . 

we get 

T = T^[T L glH Xc P c G r c P c H cp r p V R r p H pc P c G a c P c H cX g a x]- 
We now define the projection G c of Gc into c, 

G c = P c G c P c = Pc(z-H' c )- 1 P c , (B27) 



where 



P c (z -H e -E x - EpJ-^c 

H c = P C H' C P C , 
Ea = H c \g\H\c, 



(B28) 

(B29) 
(B30) 
(B31) 



Yip — H C pQ pH pc i 
so to get 

T = Tr[r L f x Hx c GlH cp r p T R g a p H pc G a c H c xg a x\- 

With a final cycle of the trace, 

T = Tr[H c xglT L ~glHx c G r c H cp g r p T R g a p H pc G a c }, 

and defining the contact resistance dressed renormalizcd 
lead injection rates T, 



T/ — H c xgx^L9xH\c 
T r = H cp g p T R g p H pc 



(B32) 
(B33) 



where the generalized Ty r contain the ballistic leads in- 
jection rates T L / R plus some other ingredients physically 
containing contact resistance. We finally get the gener- 
alized Fisher-Lee formula 

T = Tr[fiG r c f r G^}. 



5. NEGF fundamental relations 

We use Kadanoff and Baym^ notations where correla- 
tion and scattering functions G <> , are defined Her- 
mitean. Let's recall the fundamental relations 



r = i(S r - S a ) = £< - 
A = i{G r - G a ) = G< 



E>, 
G>. 



(B34) 
(B35) 



In steady-state NEGF, we have two more dynamical 
equations to get G < and G > : 



G< 
G> 



G r E < G a , 
G r E > G a , 



where £ is 



(B36) 
(B37) 



(B38) 



a sum over the leads T (T = L or R in case of two termi- 
nals) and the correlation S. E 



r/a/</> 



is the correlation 



retarded/advanced self-energy or the in/out scattering 
function. For ballistic leads T at equilibrium with their 
reservoirs one has 



T {E) = f™(E)T T (E), 

= (i-/f D (£))r T (£), 



(B39) 
(B40) 



where Tt is the injection rate of the considered ballistic 
lead T. fjP(E) is the Fermi-Dirac distribution of the 
reservoir connected to the ballistic lead under considera- 
tion. We also remind the relation 

A = G r TG a = G a TG r . (B41) 



6. Generalized Meir-Wingreen formula 

The correlation self-energy acts only in the true central 
device c, since one assumes that there are no interactions 
in the leads (ballistic and non ballistic parts). Therefore 
we have the following identity 



PcEcorvPc E corr , 



(B42) 



while, as before, for the considered terminal T = L or R; 
t = p or A we have 

P T S T Pr = S T . (B43) 

We introduce again the projectors 

P T +P C + P T , = 1, 

where P c projects on the true central device c and P T and 
P T i project on the non ballistic part of the leads. We also 
introduce the conjugated projector 

Q T = 1-P T . 

The Meir-Wingreen formula for the current irdE that 
enters in the central device through the terminal T per 
energy interval dE, is: 

ir = ^Tr[S<G>-S>Gg] 

r rvrv<^ ir v>/^ a v>/^ 11 ' v<^ a 1 

— — W\2u t Kj c U U c ~ 2^ t U c 1j Lr c J. 



%t can be divided into the coherent and non-coherent 
parts, 

%r = «T h + «T° h > (B44) 

where the non-coherent part of the Meir-Wingrcen for- 
mula is defined 

i^ coh = — Tr[£y G r c Y^ orr G a c - Y^G r c Y^ orr G a c \, (B45) 

while, using the fundamental relations of NEGF, the co- 
herent part can be written 

i-h = e Tr[r T G^r T ,G^](/| D - /£?), (B46) 

for the specific case of two terminals T and T". /|! D is 
the Fermi-Dirac distribution of lead T (always at equi- 
librium). The latter can be recognized as the ordinary 
Fisher-Lee formula when associating T = L and T' = R, 

= £ TT[ r L G r c r R G a c ](fl D - /| D ). (B47) 

To generalize the Meir-Wingreen formula we refer 
again to the the 1-c-r scheme and the associated quanti- 
ties, the Green's function of the true central device and 
the renormalized injection rates of the leads. The deriva- 
tion for the generalized Meir-Wingreen formula proceeds 
separately on the coherent and the non-coherent part. 
For the coherent part, the derivation gets back to the 
generalized Fisher-Lee formula: 

h 

For the non coherent part, we start from Eq. (|B45|) . 
Let's for the moment consider only the first term of 
Eq. (|B45|) . 

T\ = Tr^^G r c T^ orr Go\. 
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Using Eqs. (|B42l) and (|B43|) . we get 

T\ = Tr[P r E<P T G^P c £ 

corr 

PcG a c ], 

and cycling the trace and using projector properties 

T\ = r Ti[^<{P T G r c Qr}PJZ orr P c {Q T G a c P T }). 

We now use the Zanzwig-Mori relations Eqs. (|B11[) and 
(|B12[) . and the definition of the g T Green's function 
Eq. |(B25|) . so to get 

Tj = Tr[S<g;P T ff c Q T G^Q T P c S> rr P c Q T G^Q T J? c P T ^]. 

Now using the definition of H TC and H CT Eqs. (|B18|) 
and ([HT9]) . 

T\ = Ti[T,^g^H TC P c GQP c T,^ orr P c GQP c H CT g^j, 

and the definition of G c Eq. (|B27)) . 

Tx = Tv[^9rH T cGlY> orr G a c H CT g% 
we get, after cycling of the trace, at 

T x = Tr[H CT g«X<g r T H TC G r c X> rr GZ\. 

We now define the renormalized lead injection rates, 

S t = H CT g^ T ~g r T H TC , (B49) 

so that Ti reduces to 

T\ = Tr[EfGgS2 > orr G°]. 

When considering also the second term in Eq. (|B45[) . we 
get finally at our generalized Mcir-Wingrccn formula, 

■ncoh e T 1 T-rv</^' r V> C ia v>/^ r V< /° a l 

l T — — ll[2-, t U c lj corr U c lj t U c Z^ corr U c \. 
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